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Abstract 



This paper explores some previously- unrecognized consequences of Lerch's 
1905 formula for the Fermat quotient, with special attention to the sums 
s(fc, TV) — X/jl!°[tp'/;vj+i J "'^hich he introduced in this context. A general- 
ization of his result is proved in the case of composite A'^, and a new proof 
given of a sharpened result by Skula (2008) 15 when A'^ is even. We also 
sharpen the criteria given by Emma Lehmer in 1938 for a Wieferich prime 
to be simultaneously a Mirimanoff prime. 

Keywords: Fermat quotient, Wieferich prime, Mirimanoff prime 

1 Introduction 

For the Fermat quotient qp{b) = {b^^^ — l)/p (mod p) we employ wherever 
possible the briefer notation qh. Henceforth all congruences are assumed to be 
modp unless otherwise stated, and [-J signifies the greatest-integer function. 
The fact that Fermat quotients can be expressed as sums involving reciprocals 
of integers in — 1} was discovered in 1850 for the case & = 2 by Eisenstein, 
who gives q2 = + + {p-2) ~ (p^\) ' Subsequent researches 

in this direction have focused on developing equivalent results entailing fewer 
terms. Later Wieferich and Mirimanoff demonstrated the connection of Fermat 
quotients with the first case of Fermat 's Last Theorem (FLT). This connection 
retains its historical interest despite the complete proof of FLT by Wiles in 1995, 
and Skula's demonstration in 1992 that the failure of the first case of FLT 
would imply the vanishing of many similar sums but with much smaller ranges, 
which cannot be evaluated in terms of Fermat quotients. 
Lerch's remarkable formula of 1905 [T2] is 

N-qp{N)=Y,k-s{k,N), (1) 
fc=i 

where 

L JV J -1 

s{k,N)^ J2 - 

For a detailed exposition in English, see Agoh et al. ([T], pp. 32-35). In 
the case of composite N , the left-hand side of H]) is evaluated using Eisenstein's 
logarithmetic property q{ah) = q{a)+q{h). The fact that the terms in { '^^X^\ p— 

1} are the mirror-image (mod p) of those in {1, ^^""'^^ } implies that s{k, N) = 
—s{N — 1 — k,N). For ease of comparison with previous literature, in our 
final results we usually restrict k so as to be less than -^^3^, or to be of a 
particular parity; but in the proofs we use whichever form seems more intelligible 
or expressive in the given situation. 
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As Lerch himself noted, the complementarity of the terms about the middle 
of the range {l,p — 1} leads to considerable simplification of ^ above, with 
^^ (jv-i) ^j^^ vanishing for odd N, so that 

N ■ qpiN) = 

N 

- (TV - 1) • s(0, iV) - (iV - 3) • 3(1, N) - (iV - 5) • s(2, N) - . . . - s{- - 1, N) 

[Neven] 

N-qp{N) = 

-{N-l)-s{0,N)^{N-3)-s{l,N)^{N-5)-s{2,N)^...^2-s{ ^^ ^ -l,N) 

[iVodd]. (2) 

These results for the Fermat quotient are so comprehensive as to subsume all 
those previously achieved. The cases N — 2 and = 4 give, respectively, 
Glaisher's 1901 results (f7], pp. 21-22, 23) s(0,2) = -2-q2 and s(0,4) = -3 •52. 
Lerch himself {J,2\, p. 476, equations 14 and 15) pretty much explicitly writes 
out s(0, 3) = — I ■ q3 and 2 • s(0, 5) + s(l,5) = — | • 95, correcting work of 
Sylvester ([35], pp. 161-162). Although Lerch could have used his formulae to 
evaluate s(0,6), it was left to Emma Lehmer in 1938 ([IHI, PP- 356ff) to point 
out that s(0,6) = —2 • q2 ^ ^ ■ qs, and we have not found in any publication 
prior to that of Granville and Sun in 1996 ([S], p. 136) an explicit statement 
that the two instances of s{k, 12) lying at the center of the range {1, ^^""'^^ } can 
be evaluated by subtraction from known results, giving s(2, 12) = —(72 + § • 93, 
and s(3, 12) = 3 • 52 ~ f ' 93- However, it is convenient to group these instances 
among the "classical" results which were completely settled and systematized 
by Lerch's method. These are summarized in Table 1 below. 

Dilcher & Skula ([3], p. 389] report numerical investigations of all possible 
values of s(fc, N) for the two values of p for which q2 vanishes, namely 1093 and 
3511. This occurs only for the "classical" cases, eliminating the possibility that 
any other sums could be simple multiples of 52. It therefore seems likely that 
our Theorems 1 and 2 comprise essentially all the linear relations which pertain 
among sums of Lerch's type. 

2 Supplementary Notations 

Certainly, not all sums figuring in the literature of Fermat quotients can be 
reduced to Lerch's type, least of all those containing the numbers of Bernoulli, 
Euler, Fibonacci, Lucas, or Pell. However, the sums studied here are the sim- 
plest representatives of an important family of interrelated sums whose other 
members we designate as follows: 
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s'(fc, N) — terms with odd denominators in s(fc, N) 



s{N - 1- k,2N) = y s{N + k,2N) (3) 



s"{k, N) = terms with even denominators in s{k, N) = - ■ s{k, 2N) 



(4) 



s"'{k, N) = terms with denominators divisible by 3 in s{k, N) = — ■ s{k, 3N) 

3 



(5) 



s*{k, N) = s"{k, N) - s'(fc, N) = -- s{k, 2N) - - ■ s{N + k, 2N) 



(6) 



K{r, N) = terms in s(0, 1) congruent to rp mod N 



l.siN-r,N)^-j^-sir-l,N) (7) 



B{b,k,N)^ ^ -. 



(8) 



The evaluation of s'{k,N) ^ is a. simple consequence of the fact that a 
series expressed as a sum of terms of odd denominator may be condensed into 
a smaller range of unrestricted terms, as follows: 

n 2L(«-1)/2J+1 p_l (p-l)/2 

7 j j 2 i ' 

. 1 / .1 2L(p-n+l)/2J L(p-" + l)/2J ■' 

jodd jodd 

Conversely, a series expressed as a sum only of terms of even denominator may 
be simplified merely by factoring out the 2 in the denominator of its summand. 
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whence the formula for s"{k, N) The formula for s*{k, N) ^ follows imme- 
diately from these results. An important special case thereof, the simplification 
of which is effected using the Corollary below, is 

5*(0,7V) = -5(1,27V), (9) 

while s*(0, 1) = — s(l,2) = —2 ■ q2 is the classic result from 1850 of Eisenstein 
(0, p. 41), and s*(0, 2) = -s(l,4) = -q2 is in effect given by Stern ([17], P- 
185). Stern ([H], p. 184) also gives s'(0, 1) = ga, and Zhi-Hong Sun (^, p. 
288) extends this result to p^. It is easy to prove the following generalization of 
a statement in Sun and Sun ([H], p. 385): 

s{k,N) + s*{k,N) = s{k,2N). (10) 

K{r,N) ([7]), the proof of the formula for which is defered to (|15p below, 
corresponds to the K„i{s,p) of Zhi-Hong Sun [T5] and to the Kp{r, m) of Zhi-Wei 
Sun [23] ; we however omit the parameter p to simplify the notation and make it 
more uniform with the rest. Obviously k[o, 2) = s"(0, 1) and K{1, 2) = s'(0, 1). 
Zhi-Hong Sun ([ID], pp. 281, 286-288, 303) evaluates K{\,i), X(l,4), X(l,6), 
and K{—1,A) mod p"^, and K{—\,2i) mod p^, as well as some analogous sums 
over smaller ranges. 

The function s"'{k,N) ([5]) is not central to our argument, and we state 
without proof what are almost certainly the only results for it that can be 
expressed solely in terms of Fermat quotients: 

,"'(0,i)^i.,(0,3)^-i.g3 (11) 
s'" (0,2)^1.. (0,6)^-^92-^-93 (12) 

S"'{l,2) = ys{l,6) = yq2. (13) 

All of these have obvious relevance to Lehmer's problem. For an interesting 
variation mod p^ see Zhi-Hong Sun ([10], p. 288). 

The function B{b,k,N) of ([5|), while lacking a generally-accepted notation, 
has been evaluated for many values of the parameters, usually with b = ^, 2, i, 
or 3. Clearly when = 1 or — 1, it includes the family of sums here designated 
s{k,N) as a special case. It has a substantial if scattered literature beginning 
in 1912 with Bachmann [5], who gives B{^,0,1) = s(l,4) = 92, along with 
results from which the sums of the terms with odd or even denominators may 
be readily deduced. Zhi-Hong Sun (f20J, p. 311, Remark 4.1) gives in effect 
B{^,0, 1) = 4 • K{—1, 4) (mod p'^), which on Sun's own showing is equivalent to 
s(l, 4) mod p, but not mod p^. The study of the more difficult B{^, 0, 2) seems 
to have begun in 1995 with Zhi-Wei Sun [52], who evaluates it (in effect) as 
— s*(0, 2) — s*(2,4). An easier proof is given in Shan & Wang [TB], and finally 
Zhi-Hong Sun [TO] expressed it in the more concise form — s(2, 8), a simplification 
the possibility of which is implied by the definition of — s*(fc, TV) and by results 
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for s{k,N) with N — 8 discussed below. _B(2,0, 1), which is congruent mod 
but not mod p'^ to — 2s(l,4) = —2q2, is evaluated mod in Zhi-Hong Sun 
[2D] . Theorem 4.1, and, somewhat more simply, in Mestrovic [13], Theorem 1.5 
Despite the statements of both of these authors that this series is studied in 
Glaisher we canot find it mentioned therein. Z.-H. Sun ([H], pt. 2, Lemma 
2.4) gives evaluations modp, depending upon Lucas sequences, of B{b, 0, 2) and 
B{^,0,2) with b any integer. 



3 A partial generalization of Lerch's formula for 
composite A^^ 

If N is prime, then Lerch's formula ([T]) appears to be the only linear relation 
which holds among his sums of order TV. However, if N is composite, ([T| can 
be manipulated in a way which does not appear to have been previously stated. 
For example, if N is even, we have 



N ■ qp{N) ~ 2N ■ qp{^) = N ■ qp{2) 

= -{s(0, N) + s(2, N) + s(4, N) + ... + s{N- 2, N)}. 

(14) 

This coincides with a result given below, but because we want to prove a gener- 
alization of the formula (when N is composite) we cannot work directly from ([T]), 
but rather must develop a somewhat broader version of the underlying theory. 

The essential idea of Lerch's formula is that for fixed N, every pair of values 
of s{k, N) is connected with every other value through the sharing of multiples 
of each other's terms (mod p). Specifically, s{k, N) by definition contains the 
reciprocals of all j in the range { [-^^^^J , [^^J + 1}, and if M is any divisor of 
N (including N itself), and r the residue of k (mod M), then the corresponding 
values of Mj reduced modp are distributed among the sums s{rM, N) through 
s{rM + M - 1, iV) as follows: 

range of k residue of Mj (mod p) 

{0,^-1} 
{^,2.f-l} -p 
{2-^,3. 1^-1} -2p 

{{M-\)^,N-l} -{M-l)p 
That is, the values of Mj in the first row fall in — 1}, while those in the 
second row fall in {p + 1, 2j) — 1} and must be reduced by p, those in the third 
row fall in {2p + 1, 3p — 1} and must be reduced by 2p, etc. 

We can evaluate the sum of this two-dimensional array of terms in two differ- 
ent ways. First, we can collect the terms belonging to each value of Mj mod p 
as follows: 
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s(0, M) = M X {terms ee (mod M) in s(0, 1)} = M ■ K{0, M) 



s(l, M) EE Af X {terms = -p (mod M) in s(0, 1)} = M ■ K{-1, M) 



s(2, M) EE M X {terms = -2p (mod M) in s(0, 1)} = M ■ K{-2, M) 



s{M - 1, M) = M X {terms = -(M - l)p (mod M) in s(0, 1)} 

= Af •is:(-(M-l),Af). (15) 

In other words TV • K(r, N) = s(— r, A^) = — s(r ~ 1, A^), a result which coincides 
with one of Zhi-Hong Sun ( 18J, pt. 3, p. 90, Corollary 3.1). Secondly, we can 
collect the terms belonging to each value of r, as follows: 

Theorem 1 

iV TV TV 

s(0, TV) + s{—,N) + s(2 . — , TV) + . . . + si{M - l)j^,N) 

= Af{s(0, TV) + s(l, TV) + s(2, TV) + . . . + s{M - 1, TV)} 

= M ■ s(0, N/M) 



TV TV TV 

s{l, TV) + s(— + 1, TV) + s(2 • — + 1, TV) + . . . + ,s((A/ - 1)— + 1, iV) 

= M{s{M, TV) + s(Af + 1, TV) + s{M + 2, TV) + . . . + s(2A/ - 1, TV)} 

= Af-s(l,TV/Af) 
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s{—-l,N) + s{2- — -l,N) + s{3- — -l,N) + ... + s{N-l,N) 

= M{s{N - M, N) + s{N - M + 1, N) + s{N - M + 2, N) + . . . + s{N - 1, N)} 

N 

^M-s{—~l,N/M). (16) 

Now, while all the relations in (fTB]) are valid, those in the bottom half of the 
display are by symmetry simply the negatives of those in the opposite rows in 
the top half, while if M is odd the statement made in the middle row is truistic 
because each side consists of sums of terms symmetrically distributed about 
the "center line" and thus vanishes mod p. Furthermore, these relations are 
degenerate when M = N, so M must be a proper divisor of N, and thus N 
cannot be prime. 

The relationships expressed in pcp may be employed in two different ways. 
The simpler expressions at the far right, if tractable, will provide an explicit 
evaluation of the sums on the far left, although as we shall see in the cases of = 
9 and N = 18, it is not a foregone conclusion that the sums are better known for 
smaller values of TV. Alternatively, the sums on the far left and in the middle are 
homogeneous, and the congruence always admits of some simplification, both 
by direct cancelation and by use of the rule s{k, N) = —s{N — k — 1, N). 

We now examine some of the implications of (fTB)) for particular cases of M, 
reserving a discussion of particular cases of N until the end. 

3.1 The case M = 2 

The case M = 2 of (fTBl) is especially interesting. Let x — N/M; then: 

s{0,2x) + s{x,2x) = 2{s(0,2a;) + s(l,2j;)} =2-s(0,x) 

sll, 2x) + six + l,2x) = 2{s(2,2a;) + s(3,2j;)} =2-s(l,x) 

s{x - 1, 2x) + s{2x - 1, 2x) = 2{s(2x - 2, 2x) + s{2x - 1, 2x)} = 2-s{x-l, x), 

(17) 

with the usual redundancy in the second half of the display. The first row of (|17p 
yields an important corollary which will be frequently invoked in the following 
pages: 

Corollary 1 

s(0, 2x) + s{x, 2x) EE 2{s(0, 2x) + s(l, 2a;)} 

^ s(0, 2x) + 2- s(l, 2x) - s{x, 2x) = 

s(0, 2x) + 2 ■ s(l, 2x) + s{x - 1, 2x) = (18) 
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Among relations connecting the sums s{k, N) for even N which do not depend 
on knowledge of the value of some sum with smaller N, this one is noteworthy 
in that it relates only three terms, whereas Lerch's ^ involves [N/2\ terms. 
Other relations compare sums of two values of s{k,N) for even N with one or 
two values of s{k, N/2). In ^5]) . subtract both sides from 2s(0, x) + s(| — l,x). 
Then 

s(0, 2x) + s{x - 2, 2x) = 2s(0, x) + s(| - 1, x). (19) 
In (fT5|) . subtract s{0,x) from both sides. Then 



s(l, 2x) + s{x - 1, 2x) = -s(0, x). (20) 
Also, using (HI]) and ©, 

s*(0,a;) EE -s(l,2x). (21) 



Corollary 2 

Adding together the first two rows of (|17p. with the use of ^ and some manip- 
ulation we derive the relation: 



s*{0,x) + s*{l,x) = -s{l,x). (22) 

This is the explanation for the equivalence of some evaluations given by Zhi- 
Hong Sun ([H], pt. 3, Theorem 3.2, nos. 2 and 3). 

3.2 The case M = N/2 

In the first two rows of (dH), let N be even, and M = N/2. Then 



Corollary 3 



s(0, N) + s{2, N) + s(4, N) + ... + s{N ^2,N) 

N N 
= y{s(0, N) + s{l, N) + s{2, N} + ... + s{-- 1, N)} 

N 



- -5(0,2) = -iV. 92 (23a) 



s(l, N) + s(3, N) + s(5, N) + ... + s{N- 1, N) 

N N N N 

^ y ^) + + 1, ^) + + 2,N) + ... + s{N-l, N)} 

N 

= -.s{l,2) = N-q2. (23b) 
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Now let N — p ~ 1; then the sums s{k,p — 1) divide {l,p — 1} into equal pieces 
of length 1, so that: 



s'(0,l) = ---s(0,2) = g2 (24a) 

-5(1,2)^^ -5(0,2) ^-<72, (24b) 

supplying a slightly different proof of some important special cases of ^ and 
([H already stated. 

3.3 The case M = 3 

The case M = 3 of the first few rows of (fTB|) will likewise be needed below. Let 
X = N/M ; then: 

s(0, 3x) + 5(2:, 3a;) + s{2x, 3x) 

EE 3{s(0, 3a;) + s(l, 3a;) + s(2, 3a;)} 

= 3 • s(0,a;) 



1 



(p-2) - 
^,s"(0,l)^ 



s(l, 3x) + s(a; + 1, 3a;) + s{2x + 1, 3a;) 

= 3{s(3, 3a;) + s(4, 3.t) + s(5, 3a;)} 

= 3-s(l,a;). (25) 

3.4 The case M = N/3 

In the first three rows of (fTB)). let iV be divisible by 3, and M — N/3. Then 
Corollary 4 

s(0, N) + s(3, iV) + s(6, AT) + . . . + s(iV - 3, iV) 

= y {s(0, iV) + s(l, TV) + s{2, iV) + . . . + ,s(- - 1, TV)} 

N , , N , , 

= - -5(0,3) = ---g3 (26a) 



s(l,iV) + s(4,Ar) + s(7, N)- 



+ s{N-2,N) 
N 



N N N 

j{s{j,N) + s(- + 1,N) + si- + 2, N) 



+ -s(f -1,^)} 



TV 



S(l,3) EE 



(26b) 
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s(2, N) + s(5, TV) + s(8, N) + . . . + s{N - 1,N) 

N 2N 2N 2N 

^ y {^(X' ^) + + l,N) + s{— + 2,N) + ... + siN - 1, N)} 

N N 
^-.s{2,3)^--q,. (26c) 

Now let A'^ = p — 1, so that p = 1 (mod 3); then the sums s{k,p — 1) divide 
{l,p — 1} into equal pieces of length 1, so that: 

i + i + i + ...+ (^^m3)4-.s(2,3).i.,3 (27a) 
i + i + i + ...+ (^^i^(2,3)4..(l,3).0 (27b) 

i + i + i + ...+ (^^i^(0,3)4..(0,3).-i.,3. (27c) 

These formulae are an improvement upon some of Glaisher ([7], p. 18), where 
only the difference between the first and third rows is evaluated. However, 
Lerch ([12], P- 476) gives a result from which the first row is easily deduced, 
and Lehmer (|10j. p. 356) gives an equivalent result with the modulus p^. 
For p = 2 (mod 3), we obtain a parallel set of relations using ([7|): 

^ + 1 + 1 + ... + ^ K(2, 3) ^ i • s(l, 3) ^ (28a) 

i + i + i + ...+ (^^m3)4..(2,3)4.,3 (28b) 
J + i + i + ...+ (^^m3)4-.(0,3).4.,3. (28c) 

These likewise improve upon Glaisher ([7], P- 18), where only the difference 
between the second and third rows is evaluated. However, Lerch ([12], P- 476) 
gives a result from which the middle row is easily deduced. Lehmer ([10], p. 
356) gives an equivalent result with the modulus p^ . 

A few additional formulae in a similar vein but involving only terms in the 
first half of the range {l,p — 1} follow easily from the foregoing developments. 
If p = 1 (mod 6), 

\^hh---^j^2^\-<'^'^^-\-^^A-'^ (2^^) 
h\^\^---^i^2^\-'^'^'^^-\-^^~\-'^^- (2^^) 
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If p = 5 (mod 6), 



111 1 1 ^^ 2 1 , 

Y + 4 + y + + 3 -(2,6). + ,3 (30a) 

Lehmer (pXlj, p. 358) gives a result equivalent to the middle rows of these last 
groupings, with the modulus . Zhi-Hong Sun ([IS|, pt- 1, Corollary 1.12) 
gives the first two rows of each group. 

4 The remaining classical formulae 

For the sake of completing an inventory of the "classical" formulae, we state the 
following: 

If p = 1 (mod 4), 



Stern ([H], p. 187) gives the first and third rows, and Glaisher ([7], pp. 22-23) 
gives the first three rows. Lehmer (lO , p. 358) gives a result equivalent to the 
first row, mod . 
lip = 2, (mod 4), 
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Stern ([H], p. 187) gives the first and third rows, while Glaisher ([7], pp. 22-23) 
gives the first row only. Lehmer ([TU], p. 358) gives a result equivalent to the 
third row, mod . 
If p = 1 (mod 6), 



i + i + ^ + ...+ (^^m6)^i-.(5,6).i.,. + l.,3 (33a) 

^ + i + ^ + ...+ (^^i^(2,6)4..(4,6).4.,. (33b) 

i + i + ^ + ...+ (^^TO6)^i-.(3,6).i.,.-i.,3 (33c) 

i + ^ + ^ + ...+ (^^W6)4..(2,6).-i.,. + i.,3 (33d) 
i + ^ + ^ + ...+ (^^i^(5,6)4.. (1,6)4-, . 

i + ^ + ^ + ...+ (^^W6)4..(0,6).-l.,.4.,3. (33f) 

Lehmer ([10], p. 358) gives a result equivalent to the first row, mod . 
If p = 5 (mod 6), 

{4 + ^ + ...+ (4^^i^(5,6)4..(l,6)4.,. (34a) 

i4 + 3^ + ...+ (^^W6)4-.(2,6).4.,. + i.,3 (34b) 

i4 + ^ + ...+ (^^m6)4..(3,6)4.,.4.,3 (34c) 

J + ^ + ^ + ...+ (4^^i^(2,6)4..(4,6).4.,. (34d) 

i + ^ + 3^ + ...+ (^^m6)^i-.(5,6)4.,. + i.,3 (34e) 

i + ^ + ^ + ...+ (^^i^(0,6)4..(0,6).-l.,.4.,3. (34f) 
Lehmer (lOi, p. 358) gives a result equivalent to the fifth row, mod p^. 

5 Skula's sharpening of Lerch's formula when A^^ 
is even 

Up to this point, our results may be regarded as a fairly mild generalization 
of Lerch's formula. In this section, however, we give a new derivation of a 
transformation which yields significantly improved results. Consider both rows 
of ([^5a|) . and let: 
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a = terms in the first row with k < N/2 
b = terms in the first row with k > N/2 
c = terms in the second row with k < N/2 
d = terms in the second row with k > N/2. 

Now a + b, c + d are already defined in ()23ap . while by definition a + c= s(0, 2) = 
—2 ■ q2, b + d = s(l, 2) = 2 • ^2, a = —d, and c = ~b. Thus, we have enough 
information to solve a and b as follows, with a corresponding to row (a) and b 
to row (b): 



Theorem 2 



s(0, N) + s(2, N) + s(4, N) + ... + s{2- [ ^^ ^ ^^ J , N) 



<i^^ -(^.,, ,35.) 



s(l, N) + s(3, N) + s(5, N) + ... + s{2- L^^^j^J + 1, N) 

It should be noted that Skula ([IS], p. 8, Corollary 2.4) proved by a somewhat 
different technique a result equivalent to the second row; and as the sum of 
(|35ap and (|35b|) is by definition s(0, 2), the value of the first row is an obvious 
consequence of Skula's result. Nevertheless, we feel that this theorem warrants 
a closer look because much of its interest lies in the way the results generated 
by the two rows supplement one another. In the left-hand sides, the values of 
s{k,N) are simply those with k of the appropriate parity and strictly less than 
N/2. When = (mod 4), the number of terms in the left-hand sides of the 
two rows above is the same; when N = 2 (mod 4), the number of terms in the 
left-hand side of the second row is one less than that in the first row. All this 
will be clearer if we write k in its even form throughout and dovetail the values 
produced by the two rows: 





.5(2,4) 








-92 




s(0,2)=s(4,6) 






= -2 


92 




s(0,4)=s(4,8) + 


,5(6,8) 




= -3 


92 


s(0,6) 


+ s(2, 6)^.5(6, 10) 


f s(8, 10) 




= -A 


92 


s(0,8) 


+ s(2, 8)^.5(6, 12) 


f s(8,12) + 


,5(10,12) 


= -b 


92 


.s(0, 10) + ,5(2, 10) 4 


s(4, 10)=s(8, 14) 


f s(10,14) 


f ,5(12, 14) 


= -G 


92 


s(0,12) + s(2,12)4 


s(4, 12)=s(8, 16) 


f s(10,16) 


f s(12,16)-| 


-s(14, 16)^-7 


92 



(36) 
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These conditions are much stronger than those of Lerch ^ or [H)) . and are 
particularly parsimonious in the cases where N is oddly even (6, 10, 14, etc.). 
The case of TV = 6 was solved by Lehmer in 1938 [TD], and that of = 10 by 
Skula in 2008 [T5] (see below). The cases of TV = 8 and iV = 10 can be proved 
directly from the evaluations in Zhi-Hong Sun (TSl, pt. 3, Theorems 3.1 and 
3.3. respectively). 

Furthermore, if in pSap and (j35b[) we let TV = p — 1 ; then again the sums 
s{k,p — 1) divide 0,p — 1 into equal pieces of length 1, so that: 

The first row of [571 supplies a new proof of a theorem of Glaisher ([7^, p. 23, §40, 
where one of the versions of the formula is printed with a missing coefficient). 

The fact that Theorem 2 results in subsets of j of the terms in {l,p — 
1} which are evaluable in terms of Fermat quotients suggests the question of 
whether subsets of ^ or even of of the terms might be amenable to such 
treatment, as is the case for some of the individual values of s(fc, N) in Table 1. 
However, attempts to isolate a collection of non-consecutive terms s(fc, N) with 
the k in arithmetic progression and comprising only i of the terms in {l,p — 1} 
reveal that the result is so evaluable only when the values of k span the entire 
range. Unlike the formulae leading to the solutions in pSap and (|35bp . here 
the formulae relating the six pieces leave their values under-determined. Even 
should additional, undiscovered relations exist, the difficulty can be shown to 
be insurmountable in general. If we attempt to evaluate the sums of s{k,N) 
for every third value before or after the midpoint -^^^^^ we confront the sums 
s(0, 12) + s(3, 12) and s(6, 12) + s(9, 12) = -s(5, 12) - s(2, 12), and conversely 
if we attempt to evaluate the sums for every second value within a range of 
length we confront the sums s(0, 12) + s(2, 12), s(4, 12) + s(6, 12) = -2 • 

s(0, 12) - 10 • 92 - I ■ 93, and s(8, 12) + s(10, 12) = -s(3, 12) - s(l, 12). None of 
these results can be expressed solely in terms of Fermat quotients because each 
entails precisely one value of s(fc, 12) which cannot be so expressed (see Table 
2). 



6 Some consequences of our results for particu- 
lar cases of 

In what follows, we have nothing to add to the results for N = 1,2,3,4,6 (all 
included in Table 1) as surveyed by such authors as Emma Lehmer [10 and 
Dilcher & Skula Rather, we shall elaborate upon the implications of our 
theorems for certain other values of iV, especially with reference to the vanishing 
of 92 or 93, and to the interesting question, considered by Lehmer, of whether 
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they can vanish simultaneously. As is well known, the failure of the first case 
of FLT would require the exponent to satisfying both the congruence q2 = 
of Wieferich (only known solutions 1093, 3511) and the congruence 53 = of 
Mirimanoff (only known solutions 11, 1006003). Whether it is possible for the 
same number to satisfy both congruences remains an open question, but we 
have been able to sharpen Lehmer's criteria somewhat. 

First, however, we must briefly review progress in the determination of 
Lerch's sums made since his own paper of 1905. Lehmer ([TD], P- 352) made the 
perspicacious observation that for results relating to q2 which involve separating 
the terms in {l,p — 1} into ranges of 2, 3, 4, or 6 equal parts, these numbers 
"can be characterized by the fact that their totient does not exceed two." Setting 
aside the exceptional computable cases s(2, 12) and s(3, 12) previously noted, 
the significance of this distinction was revealed in 1991 when H.C. Williams 
([77]. p. 440) showed that the evaluations of s(l,5), s(l,8), and s(l,12) de- 
pend respectively upon the Lucas sequences f7p_(5^(l, —1), [/p_(2-)(2, —1), and 
C^p_( 3) (4, 1), the first and second of which correspond to the well-known Fi- 
bonacci numbers and Pell numbers; here (7) is the Legendre symbol. We shall 
say more below about the evalution of s{k, 8) and s{k, 12). Zhi-Hong Sun ((18j. 
Corollary 2.4) evaluated s*(0,9), equivalent to — s(l,18), in terms of q2 and a 
complex recurrence relation, while interestingly, the larger ranges s{k, 9) have 
still not been evaluated for any value of k. Sun & Sun ([21], p. 385) evaluate 
s(0, 10) in terms of (72 and Fibonacci numbers, and Zhi-Hong Sun presents a 
group of similar formulae involving our K{r, 10), some also dependent on q^, 
from which the values of s{k, 10) for other k can be derived ([TB], Theorem 3.1); 
also expressions for s(l,10) and s(2, 10) involving Lucas numbers (Corollary 
1.11), and some formulae involving s*{k, 15) and a variant of our K{r, 15) with 
terms of alternating sign (Theorem 3.2). 

Zhi-Wei Sun ([24 , p. 2216) points out that s{k,N) can be obtained by 
subtraction from known results for certain values of k when N — 24,40,60 
(i.e. when N has no prime-power divisors other than 2, 4, 8, 3, or 5). These 
all involve recurrence sequences such as the Fibonacci numbers, and cannot be 
evaluated solely in terms of Lerch's sums. 

6.1 = 8 

In the first row of (HH), set a; = 4. Then s(0, 8) + s(4, 8) ee 2{s(0, 8) + s(l, 8)} = 
2-s(0,4) = -6-92 which implies s(0, 8)-f 2-s(l, 8)-(-s(3, 8) = 0. In the first row of 
(dSI), set iV = 8 and M = 4, giving s(0, 8) -f s(2, 8) + s(4, 8) + s(6, 8) = 4 • s(0, 2). 
When g2 = 0, a straightforward calculation then gives 

s(0, 8) = -s(l, 8) = -s(2, 8) = s(3, 8), 

and pairwise, each of these relations is a necessary and sufficient condition for 
the vanishing of 92- While not without theoretical interest, such conditions do 
not entail fewer terms than those involving s(fc,4). 

As to the actual values of s{k, 8), Williams {^7\, p. 440) evaluates the sum 
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C/p_(|)(2,-l) = s(l,8) + s(2,8), 

where C/ is a Pell number. With the application of our Corollary 1, the right- 
hand side can be evaluated as s(0, 2) — s(0, 8) — s(3, 8) = — 2 • g2 + 2s(l,8), 
allowing the values of s(/c,8) to be obtained for all k (see Table 2). These are 
also readily obtainable from the values of K{r,8) tabulated in Zhi-Hong Sun 
Theorem 3.3). As every value of s{k,8) entails a Pell number, there is 
no reason to expect that they should be expressible as simple multiples of q2, 
and the fact that they generally cannot be is proven by the following cases 
where they vanish while q2 does not (the calculations have been extended to 
p < 1431906 without finding any further solution): 

s(0,8) p = 269, 8573, 1300709 

s(l,8) p = 29 

s(2,8) p = 193 

s(3,8) p = 23, 56993. 

6.2 = 16 

In (dlD, let M = 4. Then 

s(0, 16) + s(4, 16) + s(8, 16) + s(12, 16) = 4s(0, 4). (38) 
In the first two rows of dTT]), let x = 8. Then 

s(0, 16) + s(8, 16) = 2s(0, 16) + s(l, 16) = 2 • s(0, 8) (39a) 
s(l, 16) + s(9, 16) = 2s(2, 16) + s(3, 16) ee 2 • s(l, 8). (39b) 

The first row implies s(0, 16) + 2 • s(l, 16) - s(8, 16) = 0. When 92 = 0, 2 • 
s(0, 8) + 2 • s(l, 8) = 2 • s(0, 4) = 0, so s(0, 16) + s(8, 16) + s(l, 16) + s(9, 16) = 0. 
Adding these expressions gives 2 -3(0, 16) + 3 -5(1, 16)-|-s(9, 16) = 0, furnishing a 
criterion for the vanishing of q2 which entails only of the terms in {l,p— 1}, 
a slight improvement on that with N — 8. However, such an improvement 
does not continue for higher powers of 2, as the relationship to s(A:,4) becomes 
too tenuous. Theorem 2 generates identities in four terms which vanish when 
92 = 0. 

6.3 = 12 

We begin by noting that the conditions s(2, 12) = —q2 + 1 • (73 = and s(3, 12) = 
3 • (72 ^ § ■ 93 = 0, which Lchmcr inexplicably overlooked, provide sharp necessary 
criteria for the simultaneous vanishing of 92 and 93. Of all such criteria, these 
sums have the smallest ranges. As in the cases of the conditions s(0, 6) = 
— 2 • 92 — I • <Z3 = and s(2, 6) = — 2 • 92 + | • 93 =0, these criteria are 
certainly insufficient individually, as proven by the following cases where the 
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sums vanish although neither 92 nor 53 does (the calculations have been extended 
to p < 300000 without finding any further solution) : 



s(0,6) p = 61 

s(2,6) p = 73, 83 

s(2,12) p = 179, 619, 17807 

s(3,12) p = 250829. 



Obviously, however, since s(2,12) + s(3,12) = 2 -(72, if these two sums vanish 
together then so does 92- Because so much is already known in the case = 12, 
we shall only note further that taking M = 2 in p6p gives: 



s(0, 12) + s(6, 12) = 2{s(0, 12) + s(l, 12)} = 2 • s(0, 6) = -4 • 92 - 3 • 93 (40a) 
s(l, 12) + s(7, 12) = 2{s(2, 12) + s(3, 12)} = 2 ■ s(l, 6) = 4 ■ ga- (40b) 



Here, the more interesting relation is the second one, which gives s(l,12) = 
s(4, 12) as another necessary and sufficient condition for the vanishing of 52; 
thus (truistically) the simultaneous vanishing of q2 and would imply 



Frobenius in his paper of 1914 ([6J, p. 676) gives precisely this condition as a 
prerequisite for the failure of the first case of FLT, but his proof is not via the 
theory of the Fermat quotient. 

As to the actual values of s(fc, 12), Williams ([27 , p. 440) evaluates the sum 



where ?7 is a term in the Lucas sequence 0, 1, 4, 15, 56, .... The right-hand 
side of this is equivalent to s(0, 2) - s(0, 12) - s(5, 12) = -2 • 52 + 2 • s(l, 12), 
allowing the values of s{k, 12) to be obtained for all k (see Table 3). In addition, 
an evaluation of s(0, 12) can be recognized with some effort in Granville & Sun 
(0, p. 119). For N = 12, only s(2, 12), s(3, 12), s(8, 12), and s(9, 12) can be 
evaluated solely in terms of Fermat quotients. 

6.4 N = 24 

As recognized by Zhi-Wei Sun ([231, p. 2216), s{k, 24) can be explicitly evaluated 
by substraction from known values when k = 2, 3, 8, 9, etc., since 



s(0,12) 



s(l,12) 



s(4,12) = s(5,12). 



C/p_(|)(4, 1) = s(l, 12) + s{2, 12) + s(3, 12) + s(4, 12), 



s(2, 24) EE s(0, 8) - s(0, 12) ee s(1, 12) + s(2, 12) - s(l, 8) 
s(3, 24) = s(0, 6) - s(0, 8) = s(l, 8) - s(2, 12) 
s(8, 24) = s(2, 6) - s(3, 8) = s(2, 8) - s(3, 12) 

s(9, 24) = s(3, 8) - s(5, 12) = s(3, 12) + s(4, 12) - s(2, 8). 



(41a) 
(41b) 
(41c) 
(41d) 
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In (|17p . let X = 12. Then the relations in the third and fourth rows are: 

s(2, 24) + s(14, 24) = 2{s(4, 24) + s(5, 24)} = 2 ■ s(2, 12) (42a) 
s(3, 24) + s(15, 24) = 2{s(6, 24) + s(7, 24)} = 2 • s(3, 12), (42b) 

which will clearly vanish if it is possible for q2 and 53 to vanish simultaneously. 

6.5 = 9 

From (Eg, s(0, 9) + s(3, 9) + s(6, 9) = 3{s(0, 9) + s(l, 9) + s(2, 9)} which implies 
2 • s(0, 9) + 3 • s(l, 9) + 4 • s(2, 9) - s(3, 9) = 0, the strongest relation produced 
by Lerch's theorem other than the ones depending on the relationship with 
s(fc, 3), including those given by [25] and other cases of with n = 9, M = 3. 
Although s(fc, 9) has not been evaluated for any value of fc, it is known that 
in general it cannot be expressed as a simple multiple of 93, as proven by the 
following cases (apart from the trivial one of /c = 4) where it vanishes while (73 
does not (the calculations have been extended to p < 1043300 without finding 
any further solution): 

s(0,9) p = 677, 6691 

s(l,9) p = 151, 457, 971, 1439, 12613 

s(2,9) p = 241, 739, 37799 

s(3,9) p = 97, 58193. 

6.6 = 18 

We cannot add much to the knowledge of this little-studied case, other than to 
point out that Corollary 1 gives 

s(0, 18) + 2 • s(l, 18) + s(8, 18) = 0, (43) 
while the second row of pSj) gives 

s(l, 18) + s(7, 18) + s(13, 18) EE 6 • 92, (44) 
and the two rows of Theorem 2 give, respectively, 

s(0, 18) + s(2, 18) + s(4, 18) + s(6, 18) + s(8, 18) = -10 • 92 (45a) 
s(l, 18) + s(3, 18) + s(5, 18) + s(7, 18) ee -8 ■ 92- (45b) 

In view of the work of Zhi-Hong Sun ([IS], part 2, Corollary 2.4) in which he in 
effect evaluates s(l, 18) (apply our[^to his result), there is no reason to expect 
that s(fc, 18) would vanish with (72 or with 53. Not only does it fail to do so for 
the two known Wieferich primes and for the two known Mirimanoff primes, but 
with the possible (but unlikely) exceptions of fc = 0, 8, it cannot in general be 
expressed as a simple multiple of either of these Fermat quotients, as proven by 
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the following cases where it vanishes while neither of them does (the calculations 
have been extended to p < 3835989 without finding any further solution) : 



s(0,18) 


P = 




s(l,18) 


P = 


47, 1777, 217337 


s(2,18) 


P = 


167 


s(3,18) 


P = 


1171, 37783 


s(4,18) 


P = 


137, 251, 1087, 1301, 2111, 5749 


s(5,18) 


P = 


4177, 1581479 


s(6,18) 


P = 


108541 


s(7,18) 


P = 


149, 35267 


s(8,18) 


P ^ 





The apparent scarcity of zeroes of s(0, 18) and s(8, 18) against those for s(l, 18) 
has no obvious explanation, as all three figure in the most restrictive relation 
known (|43p , and the first two are not known to be more highly constrained than 
s(4, 18), which has more zeroes than any other value of s(fc, 18) in the range 
tested. Nor is it apparent why the distribution of zeroes in the tested range is 
so strikingly different from that of = 8. 

We believe Dilcher & Skula ( 4 , pp. 389-390) are in error when they state 
that there are zeros of s(fc, N) with p < 2000 for all values of N from 2 to 46 
except 5. 

6.7 = 5 

We shall not attempt to treat this case in any detail, for as previously noted, 
s(fc, 5) cannot generally be expressed in terms of Fermat quotients. We merely 
note that from Lerch's formula ([T]), if = 0, then 

5g5 = -4s(0,5)-2s(l,5) = 0, 

which implies that 

-2s(0,5) - s(l,5) = 0. 

From this condition and the obvious relationships between s(fc, 5) and s(fc, 10), 
it is straightforward to deduce that if (75 = 0, then 

s(0,5) EE s(4, 10) 

and 

s(l,5) = s(l,10) = -s(3,10). 

Furthermore, the results for s(fc, 10) discussed below show that if 52 = 95 = 0, 
then if either of the quantities s(0, 5) or s(l, 5) vanishes, so does the other. 
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6.8 N = 10 



An early result for this value of N was given in terms of Fibonacci numbers by 
H.C. Williams p. 369), who evaluated s*(0, 5), equivalent to — s(l,10). 

Later, the work of Zhi-Hong Sun {[TB], Corollary 1.11 and Theorem 3.1) provided 
explicit evaluations of K{r, 10), and thus indirectly of s{k, 10), for every value of 
k. Yet it is nonetheless interesting to consider the relations which pertain among 
these sums. Skula ([H], pp. 9-10) indeed made a special study of this case and 
gives s(0, 10) + 2 • s(l, 10) + s(4, 10) = 0, which corresponds to our Corollary 1 
with x = 5, and 2 • s(0, 10) + 3 • s(l, 10) + 2 • s(2, 10) + 3 • s(3, 10) + 2 • s(4, 10) = 0, 
which corresponds to our Theorem 1 with M — 5. These results can be easily 
read out of Sun's work (Theorem 3.1) although he does not state them explicitly. 

Now Zhi-Hong Sun ([TB], pt. 3, Theorem 3.2, nos. 3 and 5), proved two 
formulae which together yield the surprising result s*{l,5) = — s*(0, 3). This 
apparently anomalous relationship between sums neither of whose N values 
divides the other becomes less mysterious when rewritten in the form s(l, 10) -I- 
s(3, 10) = — s(0,6) — s(2,6), which reveals it as a direct consequence of our 
Theorem 2. When q2 = the right side of this congruence vanishes, giving 
s(l, 10) -I- s(3, 10) = as a remarkably compact condition, both necessary and 
sufficient, for the vanishing of 52- Comparison with Skula's second result then 
shows that when 52 = 0, we have also s(0, 10) -I- s(2, 10) -I- s(4, 10) = 0. If in 
addition we assume that (75 = 0, the relations simplify as follows: 

s(0, 10)= 3s(0,5) 
s(l,10)= s(l,5) 
s(2,10)=-4s(0,5) 
s(3,10)= -s(l,5) 
s(4, 10)= s(0,5) 

or, more succinctly, 

4s(0, 10) EE -6s(l, 10) = -3s(2, 10) = 6s(3, 10) = 12s(4, 10). 

The derivations arc not difficult, and we skip the details as these relations can 
be easily inferred from the explicit evaluations of K{r, 10) given in [18j . 



7 Lehmer's problem 

Although Emma Lchmer was not the first author to pose the question of whether 
q2 and 53 can vanish simultaneously, her 1938 paper remains the most impor- 
tant contribution to the subject. Indeed, there does not seem to have been 
much produced since, other that an heuristic argument against the possibility 
in Lenstra [llj . However, to the extent that Lehmer develops congruences for 
Fermat quotients to higher moduli, or derives expressions which cannot be ex- 
pressed in terms of Lerch's sums, her work is supplemented by the extensive 
writings of Zhi-Hong Sun, notably by a major recent paper on Bernoulli and 
Euler numbers |20j . 
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As previously noted, Lehmer overlooked the conditions involving s(fc, 12) 
discussed above, including the sharpest of all necessary criteria requiring only 
Fermat quotients, i.e. s(2, 12) = -~q2 + | • ga = and s(3, 12) = 3 ■ q2 ~ |g3 = 
0. We have given some comparable conditions involving N = 24. As to the 
vanishing of q2 alone, see our ([55]) . and of alone, our (|26ap : but although each 
of these results implies an infinite family of conditions, they do not appear to 
combine in any interesting way. 

8 Remark on a result of Dilcher and Skula 

Dilcher and Skula prove in [3] that the failure of the first case of Fermat's Last 
Theorem would imply 

s{k,N)=0 (46) 

for all iV < 46 and all < N. Such a failure would of course entail the vanishing 
of 92- On that assumption, adding together pSap and (|35b|) and cancelling the 
vanishing sums s(0, N) + s(l, TV), etc., we are left in the case N = 2 mod 4 with 
but a single term, s(2 • [i^^^J,7V) = s{f -1,N) = 0, and applying we 
find that all the values of s{k, N) may be expressed as multiples of values of 
s(fc, y)- In other words, when the condition is proved for an odd N and 
all k < N, the same condition is immediately proved for the case of 2A^. Thus, 
the result of Dilcher and Skula for iV < 46 automatically extends to all oddly 
even TV < 90. 

This observation complements a result of Cikanek [3], in which it is shown 
that the failure of the first case of FLT implies (|^ for all < 94 and all 
k < N. Cikanek's proof requires the additional condition (stated in §3.4 of the 
paper) that p > 5iN-if(N-2f/4^ 
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Table 1: Complete list of Lerch's sums (with k < N/2) which can be evaluated 
solely in terms of Fermat quotients 



s(0,l) 





s(0,2) 


-2-92 


s(0,3) 


-|-<Z3 


s(l,3) 





s(0,4) 


-3 • (?2 


s(l,4) 


92 


s(0,6) 


-2 • (72 - 1 • 93 


s(l,6) 


2-92 


s(2,6) 


-2 • 92 + i • 93 


s(2,12) 


-92 + 1 • 93 
3 • 92 - 1 • 93 


s(3,12) 



Table 2: Values of s(fc,8), where Un = (7p_(2)(2, —1) is a Pell number, and (7) 
is the Legendre symbol; derived from Williams [27], p. 440 

s(0,12) -4.92 - 2-{l)-f 

5(1,12) 92 + 2.(f).^ 

s(2,12) -92 + 2-il)-f 

s(3,12) 2-92 - 2.(f)-^ 



Table 3: Values of s(fc,12), where C/„ — L'p_(3j(4, 1) is a term in the Lucas 
sequence 0, 1, 4, 15, 56, . . . , and (-) is the Legendre symbol; derived from 
Williams fZ7], p. 440 



s(0,12) 


-3 


■ 92 - 


3 
2 


• 93 


- 3 


s(l,12) 




92 






+ 3 


s(2,12) 




-92 + 


3 


• 93 




s(3,12) 


3 


■ 92 - 


2 


• 93 




s(4,12) 


-3 


■ 92 






+ 3 


s(5,12) 




92 + 


3 
2 


• 93 


- 3 
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